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BACKGROUND OF THE INVENTION 
1 . Field of the Invention 

The present invention relates to a process and a device for 
displacing a moveable unit on a base. 

Said device is of the type comprising a 
controllable actuator, for example an electric motor, 
intended to give rise to a linear displacement of the 
moveable unit on the base, as well as a system which is 
formed of a plurality of elements which are brought into 
motion upon the displacement of said moveable unit. 

Within the context of the present invention, said 
system exhibits at least two different motions and comprises 
as elements which may be brought into motion, in particular: 

- said base which can be mounted elastically with respect to 
the floor, especially so as to isolate it from vibrations 
originating from said floor; and/or 

- one or more auxiliary masses, for example measurement 
supports and/or loads, which are tied elastically to the 
base; and/or 

- one or more auxiliary masses, for example likewise 
measurement supports and/or loads, which are tied 
elastically to the moveable unit. 
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When the moveable unit is set into motion, said 
elements of the system begin to move. However, especially by 
reason of the aforesaid elastic link, these elements still 
continue to move when the displacement of the moveable unit 
has terminated and when the latter comes to a stop. 

Such a continuance of the motions of said system 
is generally undesirable, since it may entail numerous 
drawbacks. In particular, it may disturb measurements, 
especially positioning measurements, which are made on the 
moveable unit or on these elements. 

Also, an object of the present invention is to 
control the moveable unit in such a way that all the moving 
elements of said system, for example the base and/or 
auxiliary masses, are stationary at the end of the 
displacement of the moveable unit. 

As regards said base, if it is mounted elastically 
with respect to the floor, it is known that, when the 
moveable unit is set into motion, during the acceleration 
and deceleration phases, it is subjected to the reaction of 
the force applied to the moveable unit by the actuator. This 
reaction load excites the base which then oscillates on its 
supports. This disturbs the relative positioning of the 
moveable unit with respect to the base, and greatly impedes 
the accuracy of the device. 
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This relative position error persists after the 
end of the displacement of the moveable unit and disappears 
only after the stabilization (which takes place much later) 
of the base. 

5 Various solutions for remedying this drawback are 

known. Some of these solutions make provision in particular: 

- to immobilize the base during the acceleration and 
deceleration phases via a disabling system, for example an 
electromagnetic disabling system, which is mounted in 

J| 10 parallel with the elastic supports. However, this known 

solution prevents the supports from isolating the base from 

'^4 the vibrations originating from the floor during said 

''^^/ 

acceleration and deceleration phases; 
M - to cancel the effect produced by the force developed by the 

W 15 actuator, by making provision for an additional actuator 

which is arranged between the base and the floor and which 

develops an additional force of the same amplitude but 

oppositely directed; or 

- to displace an additional moveable unit on the base 
2 0 according to a similar displacement, but oppositely 

directed, with respect to the displacement of the moveable 
unit, so as to cancel the inertia effects. 
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However, none of these known solutions is 
satisfactory, since their effectivenesses are restricted and 
since they all require supplementary means (disabling 
system, additional actuator, additional moveable unit) which 
increase in particular the complexity, the cost and the 
bulkiness of the device. 

Moreover, above all, these solutions implement an 
action which acts only on the base and not on the other 
elements of the system which, for their part, continue to 
move when the moveable unit is stationary. 

The object of the present invention is to remedy 
these drawbacks. It relates to a process for displacing, in 
an extremely accurate manner and at restricted cost, a 
moveable unit on a base mounted for example on the floor, 
whilst bringing all the motions to which this displacement 
gives rise to a stop at the end of the displacement, said 
moveable unit being displaced linearly according to a 
displacement which is predetermined in terms of distance and 
time, under the action of a controllable force. 

Accordingly, said process is noteworthy according 

to the invention in that: 

a) equations are defined which: 
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- illustrate a dynamic model of a system formed by 
elements, of which said moveable unit is one, which are 
brought into motion upon a displacement of said moveable 
unit; and 

- comprise at least two variables, of which the position 
of said moveable unit is one; 

b) all the variables of this system, together with said 
force, are expressed as a function of one and the same 
intermediate variable y and of a specified number of 
derivatives as a function of time of this intermediate 
variable, said force being such that, applied to said 
moveable unit, it displaces the latter according to said 
specified displacement and renders all the elements of 
said system immobile at the end of said displacement; 

c) the initial and final conditions of all said variables are 
determined; 

d) the value as a function of time of said intermediate 
variable is determined from the expressions for the 
variables defined in step b) and said initial and final 
conditions ; 

e) the value as a function of time of said force is 
calculated from the expression for the force, defined in 
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step b) and said value of the intermediate variable, 
determined in step d) ; and 
f) the value thus calculated of said force is applied to said 
moveable unit. 

Thus, the force applied to the moveable unit 
enables the latter to carry out the predetermined 
displacement envisaged, especially in terms of time and 
distance, whilst rendering the elements brought into motion 
by this displacement immobile at the end of the displacement 
so that they do not oscillate and, in particular, do not 
disturb the relative positioning between themselves and the 
moveable unit . 

It will be noted moreover that, by reason of this 
combined control of said moveable unit and of said moving 
elements, one obtains an extremely accurate displacement of 
the moveable unit in a reference frame independent of the 
base and tied for example to the floor. 

It will be noted that the implementation of the 
process in accordance with the invention is not limited to a 
displacement along a single axis, but can also be applied to 
displacements along several axes which can be regarded as 
independent , 

Advantageously, in step a) , the following 

operations are carried out: the variables of the system are 
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denoted xi, i going from 1 to p, p being an integer greater 
than or equal to 2, and the balance of the forces and of the 
moments is expressed, approximating to first order if 
necessary, in the so-called polynomial matrix form: 
A(s)X - bF 

with: 

• A(s) matrix of size p x p whose elements Aij (s) are 
polynomials of the variable s = d/dt; 

• X the vector : 

• b the vector of dimension p; and 

• F the force exerted by the motor. 

Advantageously, in step b) , the following operations are 
carried out : 

- the different variables xi of said system, i going from 1 
to p, each being required to satisfy a first expression of 
the form: 

j = r 

xi = E pi, j.y^^ , 

3 = 0 

the y^^^ being the derivatives of order j of the 
intermediate variable y, r being a predetermined integer 
and the pi, j being parameters to be determined, a second 
expression is obtained by putting y^^^=s^.y: 



xi = ( 2 pi, j.s^ )y = Pi(s) .y , 

j=0 
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- a third expression of vectorial type is defined on the 
basis of the second expressions relating to the different 
variables xi of the system: 

X = P.y 

comprising the vector P = 



- said vector P is calculated, by replacing X by the value 
P.y in the following system: 

B^.A(s) .P(s) = Op " 1 
j=p 

bp . F = S Ap, j(s) .P j(s) ,y 

3=1 

in which: 

. b"^ is the transpose of a matrix B of size px(p-l), such 
that B^b=0p-1; 

. bp is the p-th component of the vector b previously 
defined; and 

. Op-1 is a zero vector of dimension (p-1) ; 

- the values of the different parameters pi/j are deduced 
from the value thus calculated of the vector P; and 

- from these latter values are deduced the values of the 
variables xi as a function of the intermediate variable y 
and of its derivatives, on each occasion using the 
corresponding first expression. 

Thus, a fast and general method of calculation is 

obtained for calculating the relations between the variables 
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of the system and said intermediate variable, in the form of 
linear combinations of the latter and of its derivatives 
with respect to time. 

Advantageously, in step d) , a polynomial 
expression for the intermediate variable y is used to 
determine the value of the latter. 

In this case, preferably, the initial and final 
conditions of the different variables of the system, 
together with the expressions defined in step b) , are used 
to determine the parameters of this polynomial expression. 

In a first embodiment, for displacing a moveable 
unit on a base which is mounted elastically with respect to 
the floor and which may be subjected to linear and angular 
motions, advantageously, the variables of the system are the 
linear position x of the moveable unit, the linear position 
xB of the base and the angular position Gz of the base, which 
satisfy the relations: 



X = y + 
xB = 
Gz = 



(rB 


r9^ 




— + 






kBj 


V 



f mB rBre _J_^ 
kB" ^ kBke ke 



m 
kB 



ke ke 



^ rBJ mBre 
^kBke"^ kBke 



y<^' + 



mBJ 
kBke 



ke 



kB 



rB y(3) ^ y(2) 



kB 



n which: 



- m is the mass of the moveable unit; 

9 

Atty Docket No. 155603-0209 

Express Mail Label No. EL692571665US 



BJY/wrj 



- mB, kB, kO, rB, r0 are respectively the masS; the linear 
stiffness, the torsional stiffness, the linear damping and 
the torsional damping of the base; 

- J is the inertia of the base with respect to a vertical 
axis ; 

- d is the distance between the axis of translation of the 
center of mass of the moveable unit and that of the baseb- 
and 

- y*^^ ^ y*^' , y^^' and y*^' are respectively the first to fourth 
derivatives of the variable y. 

This first embodiment makes it possible to remedy 
the aforesaid drawbacks (inaccurate displacement, etc) 
related to the setting of the base into oscillation during 
the displacement of the moveable unit. 

In a second embodiment, for displacing on a base a 
moveable unit on which are elastically mounted a number p of 
auxiliary masses MAi, p being greater than or equal to 1, i 
going from 1 to p, advantageously, the variables of the 
system are the position x of the moveable unit and the 
(linear) positions zi of the p auxiliary masses MAi, which 
satisfy the relations: 
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/ n / 



X = 



TT mi 2 ri . 

n — s^+— s+1 

11 I,: 



Zl = 



i=1 V 

p r 



ki 



j=1 V '^J 



^s^ + — s + 1 



■y 



— s + 1 
ki 



■y 



in which: 

- n illustrates the product of the associated expressions; 

-mi, zi, ki and ri are respectively the mass, the position, 
the stiffness and the damping of an auxiliary mass MAi; 

- m j , kj and rj are respectively the mass, the stiffness and 
the damping of an auxiliary mass MAj ; and 



- s=d/dt. 



In a third embodiment, for displacing a moveable 



unit on a base which is mounted elastically with respect to 
the floor and on which is elastically mounted an auxiliary 
mass, advantageously, the variables of the system are the 
positions x, xB and zA respectively of the moveable unit, of 
the base and of the auxiliary mass, which satisfy the 
relations : 

X = (mAs^ + rAs + kA).{mBs^ + (rA + rB)s + (kA + kB)) - (rAs + kA)^ J.y 
xB = -My*^* 

zA = -M{rAy'^' +kAy'^>) 

in which: 

- M, mB and mA are the masses respectively of the moveable 
unit, of the base and of the auxiliary mass; 
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" rA and rB are the dampings respectively of the auxiliary 
mass and of the base; 

- kA and kB are the stiffnesses respectively of the auxiliary 
mass and of the base; and 

- s=d/dt. 

In a fourth embodiment, for displacing on a base 
mounted elastically with respect to the floor, a moveable 
unit on which is elastically mounted an auxiliary mass, 
advantageously, the variables of the system are the 
positions x, xB and zC respectively of the moveable unit, of 
the base and of the auxiliary mass, which satisfy the 
relations : 



X = [(mCs^ + rCs + kC).(mBs^ + rBs + kB)].y 
xB = [(mCs^ + rCs + kC).{Ms^ + rCs + kC) - (rCs + kC)^]y 



zC = (rCs + kC) . {mBs^ + rBs + kB) .y 

in which: 

- M, mB and mC are the masses respectively of the moveable 
unit, of the base and of the auxiliary mass; 

- rB and rC are the dampings respectively of the base and of 
the auxiliary mass; 

- kB and kC are the stiffnesses respectively of the base and 
of the auxiliary mass; and 

- s=d/dt. 
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The present invention also relates to a device of 
the type comprising: 

- a base mounted directly or indirectly on the floor; 

" a moveable unit which may be displaced linearly on said 
base; and 

- a controllable actuator able to apply a force to said 
moveable unit with a view to its displacement on said base. 

According to the invention, said device is 
noteworthy in that it furthermore comprises means, for 
example a calculator: 

- which implement steps a) to e) of the aforesaid process, so 
as to calculate a force which, applied to said moveable 
unit, makes it possible to obtain the combined effect or 
control indicated above; and 

- which determine a control command and transmit it to said 
actuator so that it applies the force thus calculated to 
said moveable unit, during a displacement. 

Thus, over and above the aforesaid advantages, the 
device in accordance with the invention does not require any 
additional mechanical means, thereby reducing its cost and 
its bulkiness and simplifying its embodiment, with respect 
to the known and aforesaid devices . 
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The figures of the appended drawing will elucidate 
the manner in which the invention may be embodied. In these 
figures^ identical references designate similar elements. 

Figures 1 and 2 respectively illustrate two 
different embodiments of the device in accordance with the 
invention . 

Figures 3 to 7 represent graphs which illustrate 
the variations over time of variables of the system, for a 
first embodiment of the device in accordance with the 
invention. 

Figures 8 to 13 represent graphs which illustrate 
the variations over time of variables of the system, for a 
second embodiment of the device in accordance with the 
invention. 

The device 1 in accordance with the invention and 
represented diagrammatically in figures 1 and 2, according 
to two different embodiments, is intended for displacing a 
moveable unit 4, for example a moveable carriage, on a base 
2, in particular a test bench. 

This device 1 can for example be applied to fast 
XY tables used in microelectronics, to machine tools, to 
conveyors, to robots, etc. 

In a known manner, said device 1 comprises, in 

addition to the base 2 and to the moveable unit 4 : 
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- supports 3, of known type, arranged between the base 2 and 
the floor S; 

- means (not represented) , for example a rail, fixed on the 
base 2 and enabling the moveable unit 4 to be displaced 
linearly on said base 2; and 

- a controllable actuator 5, preferably an electric motor, 
able to apply a force F to said moveable unit 4 with a view 
to its displacement on the base 2. 

Within the context of the present invention, the 
device 1 comprises a system SI, S2 which is formed of 
various elements specified hereinbelow and variables 
according to the embodiment contemplated, which are brought 
into motion upon the displacement of the moveable unit 4. 

According to the invention, said device 1 is 
improved in such a way as to obtain directly at the end of a 
displacement of the moveable unit 4: 

- accurate positioning of the latter in a reference frame 
(not represented) , independent of the moveable unit 4 and 
of the base 2 and tied for example to the floor; and 

- immobilization of all the moving elements of said system 
SI, S2. 

To do this, the device 1 moreover comprises, 
according to the invention, calculation means 6 which 
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calculate a particular force which is intended to be 
transmitted in the form of a control command to the actuator 
5, as illustrated by a link 7, and which is such that, 
applied to said moveable unit 4, it produces a combined 
effect (and hence combined control) : 

- on the one hand, on the moveable unit 4 so that it exactly 
carries out the envisaged displacement, especially as 
regards the prescribed duration and prescribed distance of 
displacement ; and 

- on the other hand, on said system SI, S2 so that all its 
moving elements are immobile at the end of the displacement 
of the moveable unit 4 , 

Accordingly, said calculation means 6 implement 
the process in accordance with the invention, according to 
which: 

a) equations are defined which: 

- illustrate a dynamic model of said system (for example 
SI or S2) formed by the different elements, of which 
said moveable unit 4 is one, which are brought into 
motion upon a displacement of said moveable unit 4; and 

- comprise at least three variables, of which the position 
of said moveable unit 4 is one; 
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b) all the variables of this system, together with said force 

are expressed as a function of one and the same 
intermediate variable y and of a specified number of 
derivatives as a function of time of this intermediate 
variable, said force F being required to be such that, 
applied to said moveable unit 4, it displaces the latter 
according to said specified displacement and renders all 
the elements of said system immobile at the end of said 
displacement ; 

c) the initial and final conditions of all said variables are 
determined; 

d) the value as a function of time of said intermediate 
variable is determined from the expressions for the 
variables defined in step b) and said initial and final 
conditions ; and 

e) the value of said force is calculated from the expression 
for the force, defined in step b) and said value of the 
intermediate variable, determined in step d) . 

Thus, by virtue of the invention, the force F 

applied to the moveable unit 4 enables the latter to carry 

out the predetermined displacement envisaged, especially in 

terms of time and distance, whilst rendering the elements 

(specified hereinbelow) which are brought into motion by 

this displacement immobile at the end of the displacement so 
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that they do not oscillate and, in particular, do not 
disturb the relative positioning between themselves and the 
moveable unit 4 . 

It will be noted moreover that, by reason of this 
combined effect or control of said moveable unit 4 and of 
said moving elements, one obtains an extremely accurate 
displacement of the moveable unit 4 in a reference frame 
independent of the base 2 and tied for example to the floor 
S. 

Of course, the implementation of the present 
invention is not limited to a displacement along a single 
axis, but can also be applied to displacements along several 
axes which can be regarded as independent. 

According to the invention, in step d) , a 
polynomial expression for the intermediate variable y is 
used to determine the value of the latter, and the initial 
and final conditions of the different variables of the 
system, together with the expressions defined in step b) are 
used to determine the parameters of this polynomial 
expression. 

The process in accordance with the invention will 
now be described in respect of four different systems (of 
moving elements) , 
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In a first embodiment (not represented) , the 
supports 3 are of elastic type and make it possible to 
isolate the base 2 from the vibrations originating from said 
floor S. The natural frequency of the base 2 on said elastic 
supports 3 is generally a few Hertz. Furthermore, in 
addition to the translational motion of the moveable unit 4 
controlled by the force an angular motion is created 
between the base 2 and the moveable unit 4. Specifically, in 
this case, the axis of the moveable unit 4 does not pass 
through its center of mass, the force produced by the 
actuator 5 creates a moment about the vertical axis. The 
rail is assumed to be slightly flexible and thus allows the 
moveable unit 4 small rotational motions about the vertical 
axis, which corresponds to the aforesaid relative angular 
motion between the base 2 and the moveable unit 4. 

Consequently, in this first embodiment, to 
displace the moveable unit 4 on the base 2 which is mounted 
elastically with respect to the floor and which may be 
subjected to a (relative) angular motion, the variables of 
the system are the linear position x of the moveable unit 4, 
the linear position xB of the base 2 and the angular 

position 9z of the base 2, which satisfy the relations: 
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c 



X = y + 



V 

nn 

xB = 

kB 



£8 

kB^ke 



mB rBre J ^ 



+ 



(2) 



V 



kB kBke ke ) 



+ 7^ y + 



rBJ mBre 



+ 



kBke kBke ) 



(3)_^mBJ ,4) 
kBke 



) 



ez = - d — 
ke 



J (4) ''Q (3) (2) 
— Y^^' -j Y^^' +V 

ke ke 

^mB ,4)^ XB (3)_^ (2) 
Ub kB 



in which 



- m is the mass of the moveable unit 4; 

- mB, kB, k9, rB, r9 are respectively the mass, the linear 
stiffness, the torsional stiffness, the linear damping and 
the torsional damping of the base 2; 

- J is the inertia of the base 2 with respect to a vertical 



axis ; 



- d is the distance between the axis of translation of the 
center of mass of the moveable unit 4 and that of the base 
2 ; and 



y*^' / Y^'^^ f y^^' y^""' are respectively the first to fourth 



,(3) 



.(4) 



derivatives of the variable y. 

Specifically, in this first embodiment, the 

balance of the forces and of the moments, the angle being 
approximated to first order, may be written: 



mx*^* = F 
mBxB*^' - - F - kBxB - tBxB'^' 



jez^^' =-dF-keez-ree2 



(1) 



(1) 
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It will be noted that, within the context of the 
present invention, a^^* is the derivative of order p with 
respect to time of the parameter a, regardless of a. Thus, 
for example, x^^' is the first derivative of x with respect 
to time. 

The calculation of the intermediate variable y is 

d 

achieved by putting s = — , x = P(s)y, xB = PB(s)y, 

dt 

0z = P9 (s)y and by rewriting the system (1) with this 
notation : 

ms^P{s)y = F 

(mBs^+rBs + kB)PB(s)y = -F 
(Js^+res + ke)Pe(s)y = -dF 

i.e.: 

(mBs^ + rBs + kB) PB (s) = - { Js^ + rGs + kG) PQ (s) = - ms^P(s) 

d 

and hence : 

Ps = + — s + 1 — + — s + 1 

kB kB kO ke 

PBs = — + — s + 1 

kB kO ke 

pes = -d — s^ + — s + 1 

ke kB kB 

From these expressions, we immediately deduce: 
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J 

(4)^lB (3)_^ (2) 
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(3) mBJ (4) 
kBke 



(1) 

The expression for y as a function of \ xB, 

xB'^' , Qz and ez*''* is obtained by inversion. However, this 
formula is not necessary in order to plan the trajectories 

of X, xB and Qz . Specifically, since we want a stop-stop 
displacement of the moveable unit 4 between xO at the 
instant tO and xl at the instant tl, with 

x^'ktO) = 0 = x^^*(tl) and xB(tO) = 0 = xB(tl) , xB^^^tO) = 0 = xB*'*(tl) and 

ez{tO) = 0 = ez(tl) , ez^^^(tO) = O = ez'^^(tl), with in addition 
F(tO)=0=F(tl) , 

we deduce therefrom through the aforesaid expressions (2) 



that y (to) -xO,y (tl) =xl and 



y(ti)=y^'\ti) = Y^'\t±) = Y^^\t±) = y^'^ti) = y^'^ti) = 0, 1 = 0.1 i.e. 
14 initial and final conditions. 

It is sufficient to choose y as a polynomial with 
respect to time of the form: 

y(t) = xO + (xl - xO) (a(t) )^ S a i (a(t) (3) 

i = 0 

with 
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t - to 

a(t) - 

tl - to 

and a > 7 and p > 6. The coefficients aO, ap are then 

obtained, according to standard methods, by solving a linear 
system. 

The reference trajectory sought for the 
displacement of the moveable unit 4 is then given by 
expressions (2) with y(t) given by expression (3) . 

Moreover, the force F as a function of time to be 
applied to the means 5 is obtained by integrating the value 
of y obtained via expression (3) in the expression 

F(t) = M.x^^'(t). 

In this first embodiment, we obtain: 



F(t) = M 



y + 



r rB re^ fmB rBrG J ^ 



kB^ke 



' kB kBk9 ke 



rBJ + mBrG^ (5} mBJ (6) 



V 



y' ' + y 

kBke J kBke 



(3A) 



with y(t) given by expression (3) . 

Thus, since by virtue of the device 1 the base 2 
is immobilized at the end of the displacement, it does not 
disturb the positioning of the moveable unit 4 in the 
aforesaid reference frame so that said moveable unit 4 is 
positioned in a stable manner as soon as its displacement 
ends. Moreover, since its displacement is carried out in an 
accurate manner, its positioning corresponds exactly in said 
reference frame to the sought-after positioning. 
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Represented in figures 3 to 7 are the values 
respectively of said variables y (in meters m) , x (in meters 
m) , xB (in meters m) , 9z (in radians rd) and F (in Newtons N) 

as a function of time t (in seconds s) for a particular 
exemplary embodiment, for which: 

- m=4 0 kg; 

- mB=800 kg; 

- kB=mB (5 . 27i) ^ corresponding to a natural frequency of 5 Hz; 

- rB=0. SVkBmB corresponding to a normalized damping of 0.3; 

- J=12 0 Nm corresponding to the inertia of the moveable 
unit 4 ; 

- k9 =J (10 , 271) ^ corresponding to a natural rotational 
frequency of 10 Hz; 

- rD=0.3VkeJ corresponding to a normalized rotational 
damping of 0.3; 

- d=0.01 m corresponding to the of f -centering of the moveable 
unit 4; 

- tl-tO = 0.4 s; and 

- xl-xO == 25 mm. 

The moveable unit 4 is displaced from the position 
xO at rest (xO^^^=0) at the instant tO, to the position xl at 
rest (xl^^^=0) at the instant tl. It is therefore displaced 
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over a distance of 25 mm in 0 . 4 s. To obtain this 
displacement^ as well as the immobilization (at the end of 
said displacement) of the various motions to which the 
displacement gives rise, the force F represented in figure 7 
should be applied to said moveable unit 4. This force is 
given by expression (3A) with y given by (3) for a = 7 

and P = 6. In this case, the coefficients aO up to a6 are 
given by aO=1716, al=-9009, a2=20020, a3=-24024, a4=16380, 
a5=-6006, a6=924. 

In a second embodiment represented in figure 1, 
the system SI comprises, in addition to the moveable unit 4, 
a number p of auxiliary masses MAi, p being greater than or 
equal to 1, i going from 1 to p, which are linked 
respectively by elastic links el to ep of standard type, in 
particular springs, to said moveable unit 4. In the example 
represented, p=3 . 

In this case, the variables of the system are the 
position x of the moveable unit 4 and the positions zi of 
the p auxiliary masses MAi, which satisfy the relations: 



X = 



|j — + — s + 1 



ZI = 



p f 



ki 



ki 



j=1 V N 



mj 2 n . 
— s"^ +— s + 1 



kj 



— s + 1 
ki 



(4) 
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in which: 

- n illustrates the product of the associated expressions; 
-mi, zi, ki and ri are respectively the mass, the position, 
the stiffness and the damping of an auxiliary mass MAi; 

- mj , kj and rj are respectively the mass, the stiffness and 
the damping of an auxiliary mass MAj ; and 

- s=d/dt. 

Specifically, the dynamic model of the system SI 
may be written: 

Mx'2> = F + y (ki(zi - X) + ri (zi<^' - x'^')) 
M izi*^' = ki (X - zi) + ri (x*^' - zi<^' ) , i = 1, . . . p. 

As in the foregoing, we wish to find laws of 
motion which ensure the desired displacement of the moveable 
unit 4, the auxiliary masses MAi (for example measurement 
devices and/or loads) being immobilized as soon as the 
moveable unit 4 stops. 

Accordingly, the intermediate variable y is 
calculated by the same approach as earlier and the 
trajectory of the moveable unit 4 is planned by way thereof. 

The intermediate variable y being required to 

. , d 

satisfy x=P (s) y, zi=Pi(s)y, i = l, p, with s= — , we must 

at 

have, substituting these relations into the system (5) : 
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(mis^+ris+ki) Pi = (ris+ki)P, i=l, p 

From this expression, we immediately derive 



P(s)- 



Vi=1 V i^' 



Pi = 



n 

j=1 V 



kj kj 



ri 
ki 



s + 1 



thereby proving the aforesaid formulae (4) . 

In this case, it may be demonstrated that the 
force F to be applied satisfies the relation: 



F(t) = 



i=1 



j=1 



The aforesaid formulae are verified and specified 
hereinbelow for two and three auxiliary masses MAi 
respectively . 

In the case of two auxiliary masses (p=2) , the 
model may be written: 



Mx'^* = F - k1(x - z1) - r1(x"^ - 2I"') - k2(x - z2) - r2(x*^* - z2"') 
mizf ^* - kl (x - z1) + r1 (x^''* - zf^) 
m2z2'^' = k2(x - z2) + r2(x*^* - z2*^') 



From this we immediately deduce 



X 
Zl 

z2 



.ml 2 r1 „ ,m2 2 ^^ 

— s + — s + 1 — + — s + ly 
k1 k1 k2 k2 ^ 

— s + 1 ( — s^ + — s + ly 
k1 k2 k2 

,r2 .ml 2 

— s + 1 — s'^ + — s + 1 y 
k2 k1 k1 



(6) 
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mi o ri 

i.e, putting — = Ti and — = 2DiTi, i = 1.2 : 

ki ki 

' X = y + 2{D1T1 + D2T2)y"' + (Ti^ + 12^ + 4DlD2T1T2)y'2' 

+ 2(D1T1T2^ +D2T2T1^)y'^> +(T1^T2^)y''^' 
zl = y + 2(D1 TI + D2T2)y'^' + (T2^ + 4D1D2T1T2)y'^' + (2D1TlT22)y<^' 
z2 = y + 2{D1 TI + D2T2)y*^* + (jf + 4D1D2T1T2)y*^* + (2D2T2T1^)y*^*. 

The expression for y, or more precisely the 
expressions for y, Y^^\ Y^^\ Y^^K y*^* and y^^^ / are deduced 
therefrom by inverting the system obtained on the basis of 
X, zl, z2, x^^* , zl^^^ , z2^^^ . 

We deduce therefrom that, to perform a 
displacement from xO at the instant tO to xl at the instant 
tl^ with the auxiliary masses at rest at tO and tl, it is 
sufficient to construct a reference trajectory for y with 
the initial and final conditions y(tO)=xO, y(tl)=xl and all 
the derivatives y^^^tO) =y^^^ (tl) =0 , k varying from 1 to 6 or 
more if necessary, and to deduce therefrom the reference 
trajectories of the main and auxiliary masses, as well as of 
the force F to be applied to the motor. 

In this case, the force F satisfies the relation: 



F(t) = [ ^s' + (ri + r2)s + (ki + k2)f — s' + — s + 1^ 
^ Xkl ki J 



- (ris + kl)| E|s^ + + lj - (r2s + k2)f I^s^ -h + l^ly. 



Ym2 2 r2 
— s"^ + — s + 1 
U2 k2 J 

nl 2 ri 
— s + — 
V ki ki J 
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Furthermore, the model for three auxiliary masses 



MAi (p=3) [see figure 1] , may be written, as earlier 



= F - kl(x - zl) - r1{x"' - zl*"") 



- k2(x - z2) - r2(x'^' - z2^") - k3{x - z3) - r3{x"' - z3^") 
k1(x - zl) + r1(x'^' -zl'^') 
k2(x - z2) + r2(x'^' -z2*^') 
k3(x - z3) + r3(x'^' -z3*^') 



>(i) 



(1) 



>(i)' 



mlzl'2> 
m2z2'^* 
m3z3'2' 



From this we immediately deduce: 



.ml 2 fl ,m2 2 ''2 -.,mo 2 
X = ( — s'' + — s + 1){ — s + — s + 1)( — s + — s + 1)y 



k1 ■ k2 
.m2 2 r2 



k1 

z1 = {— s + 1)(— s^+— s + 1)( 
k1 k2 k2 k3 

r2 ..,m1 2 r1 



m3 

k2 ■ k3 
m3 2 



r3_ 
k3 



^ '■3 
s + — s+1 y 
k3 

z2 = (— s + 1)(— s^ + — s + 1){— s2+— s + 1)y 
k2 k1 k1 k3 k3 

^ ,r3 „,m1 2 r1 ..,m2 2 ''2 „ 
z3 = — s + 1 — s^ + — s + 1 — s^ + — s + ly 
k3 k1 k1 k2 k2 



(7) 



We proceed as earlier in order to determine the 
values as a function of time of the different variables and 
in particular of the force F, the latter satisfying the 



expression: 



F(t) = [ (ms' + (r1 + r2 + r3)s + (k1 + k2 + k3)). 

k1 



ml 2 r1 

— s + — s + 1 

k1 



^^m2 2 r2 A(m3 2 r3 
s^ + — s + 1 s^ + 



_(ris + kl)|— s^+ — s + 1 
^ ^' k2 k2 



k2 k2 Ak3 

r2 .Vm3 2 r3 

s^ + — s + 1 



k3 



s + 1 



k3 



k3 



- (r2s + k2) 

- (r3s + k3) 



ml 2 r1 .Vm3 



■s + — s + 1 
k1 k1 

^m1 , r1 



V 



V 



s^ +— s + 1 
k1 k1 



k3 

^^m2 2 
s + 



J 

2 r3 

s + — s + 1 



k2 



k3 
r2 

k2 



s + 1 ]y. 
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Represented in figures 8 to 13 are the values 
respectively of the variables zl, z2 , 23 and F as a 

function of time t for a particular example of the 
embodiment of figure 1, zl to z3 being the displacements of 
5 the auxiliary masses MAI, MA2 and MA3 respectively. The 
variables y, x, zl, z2 and z3 are expressed in meters (m) 
and the force F in Newtons (N) . 

This example is such that: 

- M-5 kg; 

10 - ml=0 . 1 kg; 

- m2=0 . 01 kg; 

- m3 = 0 - 5 kg ; 

- kl=ml(5 .271) k2=m2 (4 .271) ^ k3=m3 (5 . 27i) corresponding to 
natural frequencies of 5, 4 and 6 Hz respectively; 

15 - rl = 0 . sVklml , r2 = 0.2Vk2m2, r3 = 0 . IsVkirrG , corresponding to 
normalized dampings of 0.3, 0.2 and 0.15 respectively; 

- tl-tO = 0.34 s; and 

- xl-xO = 4 0 mm. 

Additionally, in a third embodiment represented in 

2 0 figure 2, the system S2 comprises the moveable unit 4, the 

base 2 which is mounted elastically with respect to the 

floor S and an auxiliary mass MA which is linked by way of 

an elastic link eA of standard type to said base 2. 
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In this case, the variables of the system are the 
positions X; xB and zA of the moveable unit 4, of the base B 
and of the auxiliary mass MA, which satisfy the relations: 

X = [(mAs^ + rAs + kA).(mBs^ + (rA + rB)s + (kA + kB)) - (rAs + kA)^].y 
xB = -My*^* (8) 
zA = -M(rAy^^' +kAy*^*) 

in which: 

- M, mB and mA are the masses respectively of the moveable 
unit 4, of the base 2 and of the auxiliary mass MA; 

- rA and rB are the dampings respectively of the auxiliary 
mass MA and of the base 2; 

- kA and kB are the stiffnesses respectively of the auxiliary 
mass MA and of the base 2; and 

- s=d/dt. 

Specifically, the dynamic model of the system S2 
may be written: 

Mx^^* = F 

mBxB*'* - -F - kBxB - rBxB^'* -k{xB -zA) -rA(xB*'^ - zA^'*) (9) 

mzA^'* = kA{xB - zA) + rA(xB*^* - zA*^') 

The intermediate variable must satisfy: 

d 

X = P(s)y, xB = PB(s)y and zA = Pz(s)y with s = — . 

dt 

Substituting these expressions into (9) , we 

obtain: 
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F = Ms2p(s)y 

< (mBs^ + (rA + rB)s + (kA + kB)) PB(s) = -Ms^P(s) + (rAs + kA) Pz{s) 
(mAs^ + rAs + kA) Pz(s) = (rAs + kA) PB(s) . 

On eliminating Pz from the last equation, it 
follows that: 

[(mAs'^ + rAs + kA)(mBs^ + (rA + rB)s -^{kA + kB)) - {rAs + kA)'^ ]PB = -(mAs^ + rAs + kA)Ms 
from which we derive : 

' P = (mAs^ + rAs + kA) (mBs^ + (rA + rB)s + (kA + kB)) - (rAs + kA)^ 
^PB = -IVIs' 
Pz = -Ms^(rAs + kA) 

thus making it possible to obtain the aforesaid expressions 
(8) , 

The values as a function of time of the different 
variables, and in particular the force F, are then obtained 
as before. 

In this case, said force F satisfies the 
expression : 

F(t) = M [ (mAs^ + rAs + kA)(mBs^ + (rA + rB)s + (kA + kB))- (rAs + kA)^] y*^* . 

In a fourth and last embodiment (not represented) , 
the system is formed of the moveable unit 4, of the base 2 
and of an auxiliary mass MC which is tied elastically to 
said moveable unit 4 . 
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In this case, the variables of the system are the 
positions xB and zC respectively of the moveable unit 4, 
of the base 2 and of the auxiliary mass MC, which satisfy 
the relations: 

X = [{mCs^ + rCs + kC) .{mBs^ + rBs + kB].y 
<x3 = [ (mCs^ + rCs + kC) . (Ms^ + rCs + kC) - (rCs + kC)^ jy 
zC = (rCs 4- kC) (mBs^ + rBs + kB) .y 

in which: 

- M, mB and mC are the masses respectively of the moveable 
unit 4, of the base 2 and of the auxiliary mass MC; 

- rB and rC are the dampings respectively of the base 2 and 
of the auxiliary mass MC; 

- kB and kC are the stiffnesses respectively of the base 2 
and of the auxiliary mass MC; and 

- s=d/dt. 

Specifically, the dynamic model of this system may 
be written: 

Mx^^* - F - kC(x - zC) - rC{x*'^ - zC^'* ) 
<mBxB*'^ =-F-kBxB-rBxB^'' (10) 
mCzC*'' =kC(x-zC) + rC(x^'* -zC^'*) 
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By using, as in the foregoing, the polynomial 



representation of the variable s= — , the system (10) 

dt 



becomes : 



{Ms^ + rCs + kC)x = F + (rCs + kC)zC 

< (mBs^ + rBs + kB)xB = - F 
(mCs^ + rCs + kC)zC = (rCs + kC)x 

which, together with the expressions for each of the 
variables as a function of the intermediate variable (and of 
its derivatives) , x = P(s)y,xB = PB(s)y,zC = Pz{s)y, finally gives : 

P = (mCs^ + rCs + kCXmBs^ + rBs + kB) 

< PB = (mCs^ + rCs + kC)(Ms^ + rCs + kC) - (rCs + kC)^ 
Pz = (rCs + kC)(mBs^ + rBs + kB) 

The construction of the reference trajectories of 
y, and then of xB, zC and F is done as indicated earlier. 
In this case, the force F satisfies: 

F(t) = - (mBs^ + rBs + kB)[ (mCs^ + rCs + kc)(Ms^ + rCs + kc)- (rCs + kC)^] y. 



be described which makes it possible to determine in a 
general and fast manner the expressions defined in the 
aforesaid step b) of the process in accordance with the 
invention, for linear systems of the form: 



A method in accordance with the invention will now 



p 




(11) 
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where the Ai,j(s) are polynomials of the variable 
s, which, in the case of coupled mechanical systems, are of 
degree less than or equal to 2 and where one at least of the 
coefficients bi is non-zero. F is the control input which, 
in the above examples, is the force produced by the actuator 
5 . 

Accordingly, according to the invention, in step 
b) , the following operations are carried out: 
- the different variables xi of said system (for example SI 
or S2) , i going from 1 to p, p being an integer greater 
than or equal to 2, each being required to satisfy a first 
expression of the form: 

j = r 

xi = S pi, j.y*^^ 

j = 0 

the y*^* being the derivatives of order j of the intermediate 
variable y, r being a predetermined integer and the pi, j being 
parameters to be determined, a second expression is obtained by 
putting y^^^^s-'.y: 



XI = 



S pi, j.s^ y = Pi(s) .y 

\3 = 0 J 



- a third expression of vectorial type is defined on the basis of 
the second expressions relating to the different variables xi of 
the system: 

X=P.y 
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PI 



p 



comprising the vectors 




X 



- said vector P is calculated, replacing X by the value P.y in the 
following expressions: 

B^.A(s) .P(s) = Op - 1 



in which: 

. B*^ is the transpose of a matrix B of size px(p-l) and of rank 

p-1, such that B"^b=0p-1; 
- bp is the p-th component of the vector b; and 
. Op-1 is a zero vector of dimension (p-1) ; 

- the values of the different parameters pi/j are deduced from the 
value thus calculated of the vector P; and 

- from these latter values are deduced the values of the variables 
xi as a function of the intermediate variable y and of its 
derivatives, on each occasion using the corresponding first 
expression. 



coefficients are the polynomials Ai,j(s), i,j=l, p, i.e.: 



f3=V 

bp.F = S Ap, j(s) .P j(s) .y 



J 



The aforesaid method is now justified. 

Let us denote by A(s) the matrix of size pxp whose 
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^Al,l{s) • 


•• Al, p(s)^ 










A(s) = 










and b - 






^Ap,l(s) • 


•• Ap, p{s)^ 








,bp; 



Without loss of generality, it can be assumed that the 
rank of A(s) is equal to p (otherwise, the system is written 
together with its redundant equations and it is sufficient to 
eliminate the dependent equations) and that hp^O . There then 
exists a matrix B of size px(p-l) and of rank p-1 such that: 
B^b=0p-1 

where T represents transposition and Op-1 the vector of 
dimension p-1, all of whose components are zero. The system 
(11) premultiplied by b"^ then becomes: 

A(s)X = Op - 1 bp F = J Ap, jxj. (12) 

As indicated earlier, an intermediate variable y 
is characterized in that all the components of the vector X 
can be expressed as a function of y and of a finite number 
of its derivatives. For a controllable linear system, such 
an output always exists and the components of X can be found 
in the form of linear combinations of y and of its 
derivatives , i.e.: 
r 

xi = J^puV^* 

j=0 
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where y*'* is the derivative of order j of y with respect to 
time and where the pi,j are real numbers which are not all 
zero, or alternatively: 



xi = ( pi. js^ )y = Pi(s)y, i = I . . . , p. 

j=0 



We shall calculate the vector P{s)= 



^P1{s)^ 



vPp(s)y 



. by 



replacing X by its value P(s)y in (12) : 



B^A(s)P{s) = Op-1, bpF = X J^^^Pj^^^y- ^^^^ 

j=i 

Consequently, P belongs to the kernel of the 
matrix b'^A(s) of dimension 1, since B is of rank p-1 and A(s) 
of rank p. To calculate P, let us denote by Al(s) , 
Ap(s) the columns of the matrix A(s) and A(s) the matrix of 
size (p-l)x(p-l) defined by: 
A(s) = (A2 (s) , . . . , Ap (s) ) . 

Let us also denote by P(s) the vector of dimension 



p-1 defined by: 
^P2(s)^ 



P(s) = 



vPP(s)y 



Let us rewrite (13) in the form 
B^Al (s) PI (s) +3*^^(3) P (s) =Op-l or alternatively B^A(s)P(s)= - 
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b'^AI (s) PI (s) . Since the matrix B^A(s) is invertible, we 
have : 

P{s) - -{B'A(s))-^ B^A1{s)Pl{s) 

i.e.: 

P(s) = (co{B^A(s)))^ B^A1(s)P1{s) (14) 

det(B^A(s)) 

where co(B'^A(s)) is the matrix of the cof actors of b'^A(s) . 

From this we immediately deduce that it is sufficient 

to choose : 

P1(s) = det{B^A(s)) 

(15) 

P(s)--{co(B^A(s)))^B^ A1(s) 

this completing the calculation of the vector P(s). 

It will be observed that if the Ai,j(s) are 
polynomials of degree less than or equal to m, the degree of 
each of the components of P is less than or equal to mp. 
Specifically, in this case, the degree of the determinant 
det (B^A(s)) is less than or equal to (p-l)m and the degree 
of each of the rows of (co (b'^A (s) ) ) '^b'^ Al(s), using the fact 
that the degree of a product of polynomials is less than or 
equal to the sum of the degrees, is less than or equal to 
(p-l)m+m=pm, hence the aforesaid result. 

In all the examples presented earlier, which model 
mechanical subsystems, we have m=2 . 
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It may easily be verified that this general method 
yields the same calculations for P as in each of the 
examples already presented hereinabove. 

We shall return to certain of the examples dealt 
5 with earlier and show how the calculation of the variable y 
makes it possible to achieve passive isolation of the 
elastic modes. 

In all these examples, the trajectories are 
generated on the basis of polynomial trajectories of the 
10 intermediate value y, which are obtained through 

interpolation of the initial and final conditions. 
Furthermore, we are interested only in the particular case 
where the system is at rest at the initial and final 
instants, thereby making it possible to establish simple and 
15 standard formulae which depend only on the degree of the 
polynomial . 

In the simplest case, where the initial and final 
derivatives of y are zero up to order 4, the sought-after 
polynomial is of degree 9: 



20 



y(tO) = yO y (tO) = 0 y (tO) = 0 y (tO) = 0 y (tO) = 0 
y(t1) = y1 y (t1) = 0 y (t1) - 0 y (t1) = 0 y (t1) - 0 



which gives : 



y(t) = yO + (y1 - yO)a^ (1 26 - 420a + 540a^ - 3 1 5a^ + TOa"^ ), a = 



t-tO^ 



(16) 
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If we ask for a polynomial such that the initial 
and final derivatives are zero up to order 5, the sought- 
after polynomial is of degree 11: 

y{t) = yO + {y1 - y0)a^(462 - 1 980cj + 3465a^ - 3080a^ + 1 386a^ - 252a^) 

still with a defined as in (16) . 

If we ask for a polynomial such that the initial 

and final derivatives are zero up to order 6, the sought- 
after polynomial is of degree 13 : 

y{t) - yO + (yl - yO)a^{1 71 6 - 9009a + 20020a^ - 24024a^ + 1 6380a^ - 6006a^ + 924c^), 
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